The mode-locking transition of random lasers 
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The discovery of the spontaneous mode- locking of lasers, i.e., the synchronous oscillation of electro- 
magnetic modes in a cavity, has been a milestone of photonics allowing the realization of oscillators 
delivering ultra-short pulses. This process is so far known to occur only in standard ordered lasers 
with meter size length and only in the presence of a specific device (the saturable absorber) . Here we 
demonstrate that mode-locking can spontaneously arise also in random lasers composed by micron- 
sized laser resonances dwelling in intrinsically disordered, self- assembled clusters of nanometer-sized 
particles. Moreover by engineering a novel mode-selective pumping mechanism we show that it is 
possible to continuously drive the system from a configuration in which the various excited elec- 
tromagnetic modes oscillate in the form of several, weakly interacting, resonances to a collective 
strongly interacting regime. By realizing the smallest mode-locking device ever fabricated, we open 
the way to novel generation of miniaturized and ail-optically controlled light sources. 
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Random lasers [1] (RLs) are made by disordered highly 
scattering materials able to amplify light when exter- 
nally pumped. The simultaneous presence of structural 
disorder and nonlinearity makes these devices a fertile 
ground to connect photonics with advanced theoreti- 
cal paradigms [2] like chaos [3], non Gaussian statistics Q, 
complexity [5| and also the physics of Bose Einstein 
condensation [6]. Historically there has been a bridge 
in the RL interpretation. In pioneering experiments a 
smooth, single-peaked emission was produced by pump- 
ing finely ground laser crystals [7|, or titania particles 
dispersed in a dye-doped solution [8|, |9|. This phe- 
nomenon has been dubbed RL with incoherent feedback 
(IFRL) because it may be explained in the framework 
of the diffusion approximation [10] that neglects interfer- 
ence and treats light rays as the trajectories of random 
walking particles. However this theoretical framework 
does not explain another kind of RL exhibitting sub- 
nanometre sharp spectral peaks [lll-ll3| associated with 
high-Q resonances 1 14J - ll7| and labeled resonant feedback 
random laser (RFRL). 

Standard multimode lasers without disorder and char- 
acterized by equispaced resonances may be driven to a 
synchronous regime through the so called mode-locking 
transition[18|,ll9|, which so far has only been shown to oc- 
cur spontaneously in the presence of a saturable absorber 
and allows to generate ultra-short light pulses [20|, J21 |. 
We show that the same transition occurs in RLs and al- 
lows to lock modes of a RFRL casting its emission in the 
typical IFRL spectrum and demonstrating the inherently 
coherent nature of the random lasing phenomenon. 

The system we consider is an isolated micrometer sized 
cluster of titania nanoparticles immersed in a rhodamine 
dye solution (see supplementary information (SI)). In our 
novel setup we use the amplified spontaneous emission 
(ASE) from the surrounding dye to pump the cluster. 
The the ASE areas are defined by shaping the beam of 
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FIG. 1: The two random lasing regimes. Panel a) shows 
three normalized spectra each obtained by averaging 100 sin- 
gle shots by pumping a stripe shaped area ( 1.6 mm long). 
Top and bottom lines are retrieved for a stripe of the same 
thickness (16 \ivfi) but different orientation (15 ° tilt). Middle 
line has the same orientation of the bottom graph but dou- 
bled thickness. Panel b) shows the spectrum for a disk shaped 
pumping (diameter 1mm) for two different pump densities. 
The insets show sketches of the pumping areas. Panels c), d), 
e) and f ) show the emitted intensity distribution correspond- 
ing to the lines on the left. Pictures are retrieved by optical 
imaging of the RL emission obtained whit a pumping energy 
density of 0.1 nJ//xm 2 . The length of the white bar in figures 
e) and f) is 16/xm. 



an external (solid state) laser with the help of a reflective 



spatial light modulator. 
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FIG. 2: Cake pumping. Panel a) depicts the cluster and the 
surrounding pumped area for G = 36°. Panel b) shows S as 
a function of 0. Squares and triangles correspond to different 
pump energies for the cluster CI. Circles corresponds to the 
average of 5 measures from different clusters. The 3 dimen- 
sional graph on panel c) shows normalized spectra ( average 
over 100 shots with energy density 0.2 nJ//xm 2 ), for different 
B. Spectra are arbitrarily shifted in frequency to superimpose 
intensity maxima. AA is the wavelength shift from the most 
intense peak. 



We first consider a stripe-shaped ASE pump area of 1.6 
mm length (see figure la) and 16/im width to obtain a 
quasi-one dimensional area that acts as a strongly direc- 
tional source of ASE hitting the sample which is located 
in the centre of the stripe. The resulting spectra display 
sharp peaks. Figure la shows the average spectra over 
100 pump pulses ("shots") obtained by collecting light 
emitted off-plane from the center of the cluster. Notably 
the spectral position of the peaks remain unchanged from 
shot to shot. Dashed and continuous black lines in fig- 
ure la correspond to stripes differing by a rotation of 15° 
(see insets). Similar results are obtained for a 32/im wide 
stripe (red line in figure la). Changing the stripe orien- 
tation activates different sets of modes, as revealed by 
the change in the peaks positions. Panels c,d and e show 
the spatial intensity distribution corresponding to aver- 
aged spectra in figure la; panels lc and le correspond 
to different stripe orientations and display uncorrelated 
intensity distributions; conversely all the spots in le are 
also present in panel Id that corresponds to the same 
stripe orientation but has larger width (red and black 



continuous lines in figure la). Hence the stripe orien- 
tation affects spatial distribution of intensity and selects 
the set of activated modes. Figure lb shows the measured 
spectra when the cluster (sample "CI") is placed in the 
center of a circular pump spot (diameter 1mm) and no 
directionality is present. In this configuration spectra are 
smooth and narrow when increasing the energy, while the 
spatial intensity is homogeneously distributed (panel If). 
Having established that we can selectively excite differ- 
ent modes we proceed to study the effect of the number 
of the input k- vectors on the RL emission properties. 

To study the transition from RFRL to a IFRL we engi- 
neered a more complex pumping design (see Figure 2a). 
An area centered on the cluster is shone consisting of a 
central disk, 150 fim in diameter (to assure homogeneous 
pumping even to the largest clusters) and two symmet- 
rical wedges of much larger radius (lmm in diameter) of 
variable orientation and aperture angle (0/2) that serve 
as ASE launch pad. A single wedge configuration leads 
to the same results but proved to be hydrodynamically 
less stable. The central circle places the cluster barely 
below lasing threshold preparing the system for lasing 
once the wedge's ASE is turned on. The angular aper- 
ture O controls the number of k- vectors with which ASE 
is produced and therefore controls the number of modes 
expected to be excited. Spectra observed for small 9 
(~10°) display several very narrow ( ~ 0.05nm) peaks 
whereas large (~100°) produces a single and smooth 
RL lineshape ( ~ Anm) as shown in figure 3 below. 

To classify a RL into the IFRL or RFRL categories 
we measure its spikiness by analyzing the Fourier Trans- 
form power spectrum (FTS) of the emission spectrum. 
We define S as the high-frequency fraction of the to- 
tal FTS area, i.e. the spectral power above a frequency 
threshold, see SI). Figure 2b shows S versus at differ- 
ent pump energies for sample CI (squares and triangles) 
and averaged over 5 different clusters at high energies 
(circles). All curves display the same trend evidencing 
a transition: after a rapid growth corresponding to an 
increase in fluence and number of excited modes (ap- 
pearing on a smooth fluorescence spectrum), S reaches 
a maximum (RFRL regime) and then the spectrum be- 
comes smoother as grows, multiple peaks merge until 
an IFRL-like emission is achieved (see panel 2c). Note 
that smoothing at high is not due to averaging because 
sharp peaks are absent in single shot spectra too. 

also affects the inter-mode spectral correlation: in 
figure 3c we show intensities for a random pair of peaks 
(J(Ai = 597.2nra) and I(X 2 = 596.7nra)) for 0=18° for 
100 shots revealing an uncorrelated regime (average spec- 
tra reported in figure 3a). For 0=360° the sublte fea- 
tures present on top of the otherwise smooth spectrum 
(figure 3b) are repeatable from shot to shot (thus char- 
acteristic of the considered cluster) and the intensities 
J (A3 = 598. 4nm) and /(A4 = 598. 7nm) are strongly cor- 
related (figure 3d). Figure 3e shows the average Pearson 
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FIG. 3: The onset of a correlated random laser. Panels a) 
and b) report average spectra from cluster CI for 0=18° and 
0=360° respectively. In panel c) and d) 100 intensity pairs 
(for the chosen wavelengths see text)for 100 single shot spec- 
tra are shown for 0=18° and 6=360° respectively. Panel e) 
reports the correlation C averaged over all possible combi- 
nations of the 15 most intense peaks versus G. Error bars 
represent statistical errors from all the 105 pairs. 



correlation C (see SI) obtained from all possible pairs 
among the 15 most intense peaks (105 pairs) versus G 
for sample CI. The onset of a strongly correlated regime 
is obtained for G = 120°. The same transition has been 
observed in all the samples considered, revealing a univer- 
sal trend in which C reaches unity (within experimental 
error) when G >180°. Measurements at fixed pumping 
volume enable us to exclude the influence of the ASE 
noise (details in SI). 

In former experiments on RFRL a tightly focused 
pump spot was used to excite a limited number of modes 
thus obtaining a spectral emission displaying narrow 



spikes [15|, |22j. In our approach, for small G, we select 



modes that are strongly coupled with a directional input 
but dwell at distant positions (see figures lc, Id and le); 
in the absence of spatial overlap their mutual interac- 
tion is negligible, and the spectra obtained feature nar- 
row peaks with limited correlation (figure 3e for low G). 
Conversely, for large G, we excite a large number of spa- 
tially overlapped resonances, which result in a strongly 
correlated emission (figure 3e for large G), and a spatially 
uniform intensity distribution without hot spots (figure 
If), due to pronounced interaction between the modes. 
We reproduced these results within the framework of 



FIG. 4: Results from numerical CMT calculations. Panel a) 
and b) display the spectra obtained from 50 modes for n c =0 
and n c =10 respectively. Panel c), d) and e) represent phases 
plotted versus time for n c =3, n c =6 and n c =10. Panel e) re- 
ports the numerically calculated C (squares) and S (triangles) 
as a function of n c . 



coupled mode theory (CMT[5|, [y, LL8j), by considering a 
set of N = 50 modes at different frequencies [23], subject 
to mode repulsion [24J, [25j and excited in random initial 
conditions by an external pump pulse (see SI). In our 
model the role of G is played by the variable 2 x n c that 
is the number of resonances to which every mode couples. 

Figure 4 reports the result of our CMT calculations: 
we show average spectra in panels 4a for n c =0 present- 
ing sharp peaks and resembling an RFRL and 4b for 
n c =10 in which an IFRL-like emission including small 
features on top is retrieved. The difference between the 
two regimes becomes manifest in the time evolution of 
the modes: while the phases of the 50 weakly coupled 
(n c =3) modes (figure 4c) oscillate uncorrelated, they be- 
gin to synchronize as the coupling increases (panel 4d 
n c =6) and finally the mode-locked regime is found for 
n c =10 (figure 4e ). Note that phases are significative 
only in the time window where the pump pulse is present 
(range [0.1,0.4] in the figure, details in SI). The numer- 
ically retrieved collective parameters C and S (reported 
in figure 4f respectively as full squares and open triangles 
as a function of n c ), agree with the experimental results. 

In conclusion, we developed a new mode-locking tech- 
nique based on a pumping scheme that enables to select 
the number of activated modes in a random laser. We 



demonstrate that RLs may be prepared in two distinct 
regimes, which we can turn on by controlling the shape 
of the pumping area. When pumping is nearly unidirec- 
tional, few (low interacting) modes are turned on and 
appear as sharp, uncorrelated peaks in the spectrum. By 
increasing the number of spatial directions hitting the 
lasing cluster, many resonances intervene generating a 
smooth emission spectrum with a high degree of correla- 
tion, corresponding to the onset of a mode-locked collec- 
tive state. Only by assuming a transition from a set of 
nearly independent resonances to a regime with strongly 
interacting and mode-locked resonances, it is possible to 
explain the passage from uncorrelated spectral peaks to 
strongly correlated collective line narrowing. This exper- 
iment not only unveils the intimate and unique nature 
of random lasers and the universality of mode locking 
process paving the way for a new generation of miniatur- 
ized optical devices with engineered and tunable spectral 
emission, but also lays the foundations for a bridge be- 
tween disordered photonics and the statistical physics of 
complex systems. 
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